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We have measured the Cooper Pair Transistor (CPT) in a tunable electromagnetic environment 
consisting of four one-dimensional SQUID arrays. The transport properties of the CPT in the high 
impedance limit, Z^nv S> Rq ~ 6.45 kf2, are studied for different ratios of the Josephson coupling 
energy to the charging energy. As the impedance of the environment is increased, the current- 
voltage characteristic (IVC) of the CPT develops a Coulomb blockade of Cooper pair tunneling and 
the measured IVCs agree qualitatively with a theory based on quasicharge dynamics for a CPT. 
Increasing the impedance of the environment induces a transition in the modulation of the IVC with 
the gate charge from e-periodic to 2e-periodic. 

PACS numbers: 73.23.Hk, 73.40.Gk, 74.50.-|-r 



I. INTRODUCTION 

The interplay between charge quantization and su- 
perconductivity offers the possibility to study macro- 
scopic quantum phenomena in small-capacitance Joseph- 
son junction circuits. One of the most extensively stud- 
ied devices in this context is the Cooper Pair Transistor 
(CPT), which consists of two Josephson junctions in se- 
ries, with a gate capacitively coupled to the small island 
formed between the two junctions. Such single island cir- 
cuits are fundamental building blocks for quantum elec- 
tronic circuits based on the number-phase uncertainty of 
the Cooper pair condensate, and have been the object 
of growing interest due to their potential to realize new 
types of quantum limited measurements and solid state 
quantum bits. Essential for all these applications is the 
ability to inhibit the excitation of quasiparticles which 
give rise to dissipation (noise), and thereby decoherence. 
For the CPT this means that the charge parity of the 
island must be well controlled, so that quasiparticles, or 
single unpaired electrons, are not transferred to or from 
the CPT island. Here we study the CPT in a tunable 
high frequency electrodynamic environment, where the 
DC current-voltage characteristic (IVC) of the CPT ex- 
hibits a Coulomb blockade of Cooper pair tunneling. By 
tuning the impedance of the environment, we demon- 
strate how the parity of the island is controlled by the 
environment. 

CPTs have been most extensively studied in a high fre- 
quency environment with low impedance, ^ Rq = 
h/Ae^ ~ 6.45 kSl, where the gate voltage modulates the 
critical current flowing through the device. At low tem- 
perature and in the absence of quasiparticle excitations, 
this modulation should be periodic in gate voltage with 
period 2e/Cg, where Cg is the gate capacitance. 2e- 
periodic modulation is a clear sign that charge trans- 



port through the device is due to Cooper pairs, which do 
not change the parity of the island. However, in many 
experiments this parity effect is spoiled by the trans- 
fer of quasiparticles to or from the island, resulting in 
an e-periodic gate modulation, or modified 2e-periodic 
modulation, as the unpaired electron occupies the island. 
For applications like the coherent control of a quantum 
bit, these quasiparticles have to be avoided, and sev- 
eral experimental studies have been carried out in or- 
der to understand and control this parity effect. Normal 
metal leads, or small normal metal "traps", contacting 
the superconducting source and drain close to the tunnel 
junctions are thought to help bring long-lifetime out-of- 
equilibrium quasiparticles into thermal equilibrium, and 
thereby enhance the parity effect P, Q- However, in 
many experiments where no traps were used, a parity 
effect was demonstrated H, ^ IE 111 ■ Aumentado et al. 
showed that the lifetime of the unpaired quasiparticle on 
the island could be dramatically reduced by fabricating 
the island of the CPT with a larger superconducting en- 
ergy gap than that of the source-drain. 

Here we present an experimental study of the CPT in 
a high impedance environment, 2^cnv 3> Rq, where far 
less experimental work has been done due to the diffi- 
culty in constructing such an environment. In a high 
impedance environment, the IVC of the CPT shows not 
a supercurrent, but rather a Coulomb blockade of Cooper 
pair tunneling. The Coulomb blockade is modulated 
with the gate voltage and, in the absence of quasiparti- 
cles, this modulation is 2e-periodic. Previously, the high 
impedance environment has been achieved with small on- 
chip resistors {R ^ Rq) located close to the CPT, which 
have enabled the observation of a Coulomb blockade of 
Cooper pair tunneling with 2e-periodic dependence on 
the gate charge 0, U More recently, one-dimensional 
SQUID arrays have been used to bias a CPT which is 
the technique used in the present study. The great ad- 
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FIG. 1; SEM micrograph of a CPT (in the center) biased 
by four SQUID arrays. Two SQUID arrays are employed 
to apply a symmetric bias, while the other pair is used to 
measure the voltage across the CPT. 



vantage of the SQUID array bias is that the impedance of 
the environment can be tuned in situ, without affecting 
the CPT. In this way the role of fluctuations due to the 
environment can be easily identified. 

In this article we study the IVC of the CPT as the 
environment impedance is changed over many orders of 
magnitude. Starting from a superconducting-like IVC 
at low impedance, a well defined Coulomb blockade of 
Cooper pair tunneling with a back-bending IVC devel- 
ops as the impedance of the environment is increased. 
This back-bending IVC is due to the coherent tunnel- 
ing of single Cooper pairs and is a clear indication that 
Zenv ^ Rq- We show that the gate voltage depen- 
dence of the IVC changes from e-periodic to 2e-periodic 
as the environment impedance is increased. The high 
impedance environment suppresses quasiparticle tunnel- 
ing rates, thereby restoring the even parity of the island. 
These observations indicate new approaches for control- 
ling out-of-equilibrium quasiparticles in superconducting 
quantum circuits. 



II. FABRICATION AND MEASUREMENT 

Figure ^ shows a SEM micrograph of a CPT sample 
biased by four SQUID arrays. The Al/AlOx/Al tunnel 
junctions are fabricated on a Si02/Si substrate using e- 
beam lithography and the standard shadow evaporation 
technique Four arrays, each consisting of 70 SQUIDs 
in series, connect the CPT to the bias and measurement 
circuitry. The SQUIDs are designed to have two identical 
junctions in order to suppress the Josephson coupling in 
the array as low as possible. The nominal area of each 



junction in a SQUID is 0.03 /im^, and the loop area is 
0.18 /Ltm^. Each CPT junction has a nominal area of 
0.01-0.02 yLtm-^. To provide an effective high impedance 
environment, the arrays have to be located close to the 
CPT, so as to avoid the shunting capacitance of the leads. 

The SQUID arrays have high impedance due to the 
Josephson inductance of the SQUID Lj — h/2elc, which 
becomes infinite as the critical current of the SQUID, 
/c, is suppressed to zero with an external magnetic flux. 
This inductance is also non-linear, as it depends on the 
phase across the SQUID, and therefore a linear circuit 
model can not give a complete description of the electro- 
dynamics of these arrays. However, for current much less 
than the critical current, a linear approximation can be 
made, which shows that the infinite array behaves like 
a transmission line with a gap, or region of frequency 
where no propagating modes exist. A finite array will 
have a set of discrete modes, and for arrays such as the 
ones used in this study, the lowest frequency mode has 
/ ~ 30 GHz US ini. The impedance of arrays such 
as these has not yet been systematically characterized 
due to the very high frequencies involved. However, we 
have shown that these arrays do provide an environment 
which induces a Coulomb blockade of Cooper pair tun- 
neling and Bloch oscillations in a single Josephson junc- 
tion As the magnetic flux through the SQUID 
loops is increased, the arrays undergo a quantum phase 
transition 0, ^3 where the measured zero bias resis- 
tance of the array increases several orders of magnitude 
(50 k^l < Rq < 1 Gfl). A key aspect of this type of 
experiments is that the arrays can be tuned in this tran- 
sition region to a point where they provide a high enough 
impedance at high frequencies, while not having too large 
Rq. If Rq becomes too large, the input offset current of 
the voltage amplifier will inhibit the measurement of the 
differential voltage across the CPT. 

The chip with the arrays and CPT is bonded to a 
printed circuit board with connectors, which is mounted 
in a RF-tight box on a dilution refrigerator with a base 
temperature of 15 mK. The cryostat leads are twisted 
constatan pairs, which act as rather good microwave fil- 
ters attenuating 35 dB/GHz up to 2 GHz, above which 
a background level of -80 dB is measured with an open 
cryostat. Additional room temperature RC filters were 
used to suppress lower frequency noise when necessary. 
No microwave filters between the sample and twisted 
pairs were implemented in these measurements. How- 
ever, in the high impedance case, the arrays themselves 
are good filters, protecting the CPT from electromagnetic 
fluctuations coming from the bias side of the arrays. 

The IVC of the CPT is measured in a four point config- 
uration with the CPT symmetrically biased through one 
pair of SQUID arrays, while the other pair is used to mea- 
sure the voltage across the sample. The current is mea- 
sured with a trans-impedance amplifier (modified Stan- 
ford Research System 570) on one of the biasing leads. 
The voltage across the CPT is measured with a high in- 
put impedance differential voltage amplifier (Burr-Brown 
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FIG. 2: Schematic of a Cooper pair transistor. The en- 
vironment is represented by the impedance Zcnv When 
Zenv Rq {Zcnv 3> Rq) the CPT is embedded in a low 
(high) impedance environment. Qi(2) indicates the charge 
and (^1(2) the phase difference across each junction. Ii(2) rep- 
resents the current flowing through junction 1(2) and Ig the 
displacement current from the island to the gate electrode. 



INAllO). The bias voltages ±14 are applied symmetri- 
cally with respect to ground in order to avoid a common 
mode drain-source voltage. In this way, a polarization 
charge can be induced on the island only by the voltage 
applied to the gate and not from a bias-induced common 
mode voltage ac ting across the stray capacitance of the 
CPT to ground [ij 

The parameters of the measured samples are listed in 
Table 1. The normal state resistance per junction Rn in 
the CPT is calculated from the slope of the IVC at large 
bias voltage, assuming that both junctions are identical. 
The Josephson energy is Ej = {Rq/Rn){^/2), where 
A « 200 ^eV is the superconducting energy gap of Al. 
The capacitance of one junction in the CPT is calculated 
from the junction area as measured from SEM pictures, 
using the specific capacitance 45 fF //im^ . From the ca- 
pacitance C, the charging energy of one junction in the 
CPT is calculated as Ec = e^/2C. For all the measured 
samples, Cg ~ 10 aF ^ C, so that the total island ca- 
pacitance is given by Cs = Ci + C2 + Cg « 2C. In 
table 1, Vt indicates the measured maximum threshold 
voltage for Cooper pairs, as determined from the IVC of 
the CPT when Z^nv 3> Rq. V^ap refers to the maximum 
value of the measured Coulomb gap for single electrons, 
estimated from the offset of the normal state IVC, taken 
in high magnetic fields where superconductivity is com- 
pletely suppressed. 



III. CURRENT- VOLTAGE CHARACTERISTIC 
OF THE COOPER PAIR TRANSISTOR IN A 
HIGH IMPEDANCE ENVIRONMENT 

Most experimental and theoretical studies to date 
have considered the case of a CPT embedded in a low 
impedance environment, Zcm Rq [H 13. 13. 13. IHIR IT6| . 
A few studies have been made for the case of a high 
impedance environment, Z^nv ^ Rq Q,H,I1|' Below we 
review the theoretical description of the CPT for the low 
impedance case, showing how it can be transformed to 



a description appropriate for the high impedance case. 
We then give a qualitative description of the IVC and 
Bloch oscillations for the later case, comparing this with 
experimental results. 

We consider the circuit depicted in Fig. [3 consisting 
of a symmetric CPT with Eji = Ej2 = Ej and Ci = 
C2 — C. The charge induced by the gate on the island is 
Qg — CgVg, whcrc Vg IS thc voltage applied to the gate 
electrode. Qi(2) indicates the charge and (/)i(2) the phase 
difference across each junction. Let cf) = (f>i + (f)2 he the 
bias phase across the CPT and (pg — {(j>i — <p2)/'2 the 
phase of the island. The phases (j) and (f)g are conjugate 
respectively to the charge transferred across the device 
Q<P — {Qi + Q2)/2 and the total island charge Q = Qi~ 
Q2 + Qg and obey the commutation relations [Q^, (j)] — 
[g, cbg] = ~2ie. 

For the moment we neglect the contribution of quasi- 
particles by assuming that the superconducting energy 
gap A is large compared to all relevant energy scales. 
The Hamiltonian of a CPT with Cg ^ C biased by a 
current lb = Vb/Zcnv can be written as follows 



H = 8Ec 



Q<t> 

2e 



+ 2Ec 



Q + Qg 



2e 



-2Ejcos{(p/2)cos{(pg) 



2e 



(1) 



The first term of the Hamiltonian corresponds to the 
charging energy associated with the total charge across 
the series combination of the two junctions of the CPT. 
The second term is the electrostatic energy associated 
with the island which sees the two junctions in parallel. 
The third term is the effective Josephson energy of the 
CPT, which, due to its dependence on (pg, can be mod- 
ulated by the gate voltage. The last term represents the 
coupling between the bias current 7^ and the phase (j). 

In a low impedance environment, Zcnv ^ Rq, quan- 
tum fluctuations of the phase are suppressed and the bias 
phase (f> behaves as a classical variable. Due to the low 
impedance. Cooper pair tunneling at zero bias voltage 
dominates the IVC of the CPT which shows the typi- 
cal superconducting behavior characterized by a super- 
current branch at vanishing voltages. When Ej < Ec, 
charging effects are observable and result in a 2e-periodic 
modulation of the critical current with the gate charge. 
Thus, when Zenv Rq, the CPT behaves as a single 
Josephson junction with a gate charge dependent critical 
current 

Let us turn to the case of a high impedance environ- 
ment, Zcnv ^ Rq, where quantum fluctuations of the 
phase are strong and the conjugate charge now becomes 
a classical variable. For this case, a quantitative theory 
of the single Josephson junction in terms of Bloch energy 
bands has been developed [TtL ITsI I . The dynamics of the 
junction can be described by the quasicharge q — J Idt, 
where / is the current flowing through the junction. The 
quasicharge evolves in time according to the Langevin 
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TABLE I: List of a few characteristic parameters of the measured CPT samples. 





Rn (kO.) 


C (fF) 


Ej ifieV) 


Ec (/ieV) 


Ej/Ec 


Vt ifiY) 




CPT 6 


5.2 


0.75 


125 


106 


1.2 


15 


35 


CPT 7 


21 


0.6 


30 


133 


0.22 


40 


110 


CPT 8 


1.5 


0.8 


430 


100 


4.3 


5 


/ 


CPT 9a 


11.5 


0.55 


56 


145 


0.38 


34 


60 


CPT 9b 


8.2 


0.55 


80 


145 


0.55 


26 


45 




FIG. 3: Lowest energy Bloch band for a CPT with Ej /Ec < 
1. For fixed gate charges, the quasicharges qi and q2 follow 
a diagonal trajectory. The solid and dotted lines indicate 
respectively the quasicharges trajectories for Qg = 2ne and 
Qg = (2n + l)e. 



equation 



(2) 



where V{q) is the voltage across the junction and SI 
indicates a random noise component of the bias cur- 
rent, induced by the real impedance of the environment 
R — Re[Zenv] ^ Rq at finite temperature. The volt- 
age V{q) is given by the derivative of the lowest Bloch 
energy band, V{q) = dEo/dq. The IVC of the junction 
can be calculated analytically from the Langevin equa- 
tion 0, 0, 01 and shows a voltage peak near zero cur- 
rent corresponding to the Coulomb blockade of Cooper 
pair tunneling, and a back-bending region at finite cur- 
rents. This region of negative differential resistance is 
the indication of the Bloch oscillations, which occur with 
frequency fs = (/)/2e. 

The natural extension of the Bloch band picture from a 
single Josephson junction to a CPT leads to a description 
of the quasicharge dynamics in two-dimensional Bloch 
energy bands. In Fig.|21the lowest energy band of a CPT 
with Ej I Eq <C 1 is plotted in a space spanned by the two 
quasicharges qi and (72, which describe the dynamics of 



the two junctions of the CPT. This description becomes 
clear if we rewrite the Hamiltonian (^) in terms of the 



variables Qi(2) and ' 
reads 



h(2) 



For h = the Hamiltonian 



H = AE, 



c 



2e 



AEr 



W2 

2e 



- £:jcos((/)i) - £:jcos((/)2). 



Q2)Qg 

(3) 



This Hamiltonian is that of two single Josephson junc- 
tions, iJi(2) = 4£;c(Qi(2)/2e)^ - Ejcos{(j)i(2)), coupled 
by the term containing the gate charge Qg. For small 
coupling, Cg ^ C, the behavior of a CPT in a high 
impedance environment is described by the quasicharges 
<Zi(2) = / Ii(2)dt, where /i(2) indicates the current flow- 
ing through junction 1(2) (Fig.|21). In this case the lowest 
energy band of the CPT is approximately given by the 
sum of the lowest energy band of the two single junctions, 
that is £^0(91, 92) ^ EQ{qi) + Eo{q2). 

The conservation of current results in the constitu- 
tive relation Ii — I2 + Ig, where Ig indicates the dis- 
placement current flowing from the island of the CPT 
to the gate. From this relation follows the condition 
qi ~ q2 — Qg, which, for a fixed gate charge, constrains 
the quasicharges to a diagonal trajectory in the two- 
dimensional Bloch band (Fig. UJ. A change in the gate 
voltage induces a shift in the trajectory of the two qua- 
sicharges, as shown in Fig. O In analogy with the one 
dimensional case, the voltage across the CPT can be 
expressed via the derivative of the lowest energy band 
along the diagonal direction as V{q) — dEo(qi,q2)/dq, 
where q = {qi + 92)/2 defines the quasicharge of the 
CPT. The maximum Coulomb blockade voltage of the 
CPT, or critical voltage Vc = max[di?o('?ij 92)/^?], is a 
2e-periodic function of the gate charge. The dependence 
of the critical voltage on the gate charge is determined 
by the slope of the energy band according to the Bloch 
band theory @. The critical voltage is maximum when 
Qg = 2ne, where the slope of the lowest energy band is 
maximized (solid line in Fig. ^ , while the minimum crit- 
ical voltage is obtained for Qg — {2n + l)e (dotted line 
in Fig. El . 

Thus, when Zcnv ^ Rq the CPT behaves as a sin- 
gle Josephson junction with a gate charge dependent 



critical voltage. In the case Qg 



2ne, the energy- 



charge relation of the CPT is identical to that of a sin- 
gle Josephson junction and the voltage V{q) across the 
CPT has a 2e-periodic dependence on the quasicharge. 
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From the band diagram one can see that as the gate 
charge is shifted, the dependence of the voltage on the 
quasicharge is modified. When Qg = {2n + l)e, the 
Bloch oscillations in each junction arc out of phase, re- 
sulting in an oscillation of the voltage across the CPT 
with the quasicharge having frequency / = e/(/) 8]. For 
Ej > Ec, the lowest Bloch band becomes more sinu- 
soidal along the quasicharge axes, and the diagonal tra- 
jectory at Qg = (2n + l)e becomes flat for a symmetric 
CPT. In this case the Coulomb blockade for Cooper pairs 
can be modulated to zero and the out of phase Bloch os- 
cillations cancel one another. 

The IVC of a CPT in the high impedance environment 
shows the same features as that for a single Josephson 
junction and can be calculated from the Langevin equa- 
tion for the quasicharge ||2J). We can picture how the 
Bloch oscillations occur in the CPT, or single Josephson 
junction with tunable critical voltage, by simply consid- 
ering a non-linear capacitor biased by a voltage source 
Vb through a resistor R 3> Rq ■ In Fig. 0] the generic 
shape of the IVC is shown, with the time evolution of 
the junction voltage sketched for four different operating 
points of the circuit. We consider the case Ej < Eq, 
where the non-linear capacitor has a periodic, saw-tooth 
form of the voltage-charge relation, given by the deriva- 
tive of the lowest energy band along a diagonal trajectory 
with a maximum critical voltage. When the bias voltage 
is below the critical voltage, Vb < Vc (operating point 
a), the Langevin equation ^ will have a stationary so- 
lution with / = dq/dt — 0, and the non- liner capacitor 
will simply charge up to the value Vb with no voltage 
oscillation. As the bias is increased so that Vf, > Vc, 
time dependent solutions to the Langevin equation are 
possible, which describe overdamped, forced Bloch oscil- 
lations in the junction. These oscillations result in a peri- 
odic charging and discharging of the non-linear capacitor, 
which lead to a time average voltage {V) < Vc (operat- 
ing point b). As the bias voltage is further increased, the 
Bloch oscillations increase in frequency as higher current 
is forced through the junction. With faster Bloch oscil- 
lations, the finite charging time due to the bias resistor 
becomes negligible, and the voltage oscillations become 
more symmetric around zero, with {V) (operat- 
ing point c). Increasing the bias voltage further would 
cause back-bending of the IVC to dissipation-less trans- 
port, or DC current with zero average voltage. However, 
with higher frequency Bloch oscillations, Zener tunneling 
to higher energy bands takes place, resulting in a tem- 
porary build-up of the junction voltage to a level close 
to the superconducting energy gap 20]. In this case, 
quasiparticle tunneling can no longer be neglected, and 
Josephson-quasiparticle tunneling cycles, which are dissi- 
pative processes, can occur. As the Zener tunneling prob- 
ability increases with the bias current, the time average 
voltage across the junction increases (operating point d). 
The magnitude of the Zener current depends on the en- 
ergy gap between the lowest and first Bloch energy band, 
which in turn depends on the ratio Ej / Ec 113, 0| . 




<V> (e/C) 



v(t), 



Vc 
Vbl- 





FIG. 4: A pictorial description of the Bloch oscillations and 
back-bending IVC of a CPT, or single Josephson junction 
with tunable critical voltage Vc, current biased by a high 
impedance resistor. Iz indicates the Zener current. The time 
dependence of the voltage across the junction is shown for 4 
different operating points in the IVC. The shape of the Bloch 
oscillations refers to the case of a junction with Ej <C Ec- 



We find qualitative agreement with this theoretical pic- 
ture in our experimental IVCs of CPTs biased by SQUID 
arrays. Figure shows the IVC of two biasing SQUID 
of sample CPT 9a measured in series at three values of 
the magnetic field, between zero and $o/2. Here we see 
that the arrays change their zero bias resistance Rq over 
several orders of magnitude as the Josephson coupling of 
the SQUIDs is suppressed. In Fig. 03 the correspond- 
ing IVCs of the CPT are plotted for the same magnetic 
fields. At low Ro of the arrays, the IVC of the CPT 
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FIG. 5: (A): IVC of the two biasing SQUID arrays of sample 
CPT 9a for different magnetic fields. The measured values for 
_Ro are; Light-gray curve 7?o ~ 50 k57, dark-gray curve _Ro ~ 
500 kSl, black curve i?o « 20 Mfi. (B) Corresponding IVCs 
of CPT 9a. 



shows a superconducting-like feature that disappears as 
the Coulomb blockade region becomes progressively more 
defined for increasing values of i?o- When i?o ~ 20 Mil, 
the Coulomb blockade for Cooper pairs is fully developed 
and is followed by a region of negative differential resis- 
tance. This back-bending IVC is the hall mark feature of 
a Coulomb blockade of Cooper pair tunneling, and clear 
evidence that a high impedance of the environment has 
been achieved. Figure shows the IVC of three different 
CPT samples with different values of E,j/Ec- In each 
case the SQUID arrays are tuned to approximately the 
same value of i?o where the back-bending IVC is well 
defined. Figure El clearly demonstrates qualitative agree- 
ment with the theory described above. With increasing 
Ej/Ec ratio, the critical voltage decreases, as the lowest 
Bloch band develops more shallow minima. We also see 
that the Zener current, becomes larger with increasing 
Ej/Ec ratio, as the gap to higher energy bands becomes 
larger. 
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FIG. 6: The IVC for three different CPTs having approxi- 
mately the same value of the zero bias resistance of the SQUID 
arrays {Ro ~ 20 MQ). Black curve: CPT 9a, Ej/Ec = 0.38, 
Light-gray curve: CPT 6, Ej/Ec = 1.2 Dark-gray curve: 
CPT 8, E.,/Ec = 4.3. 



IV. GATE VOLTAGE MODULATION AND 
IMPEDANCE DEPENDENT PARITY EFFECT 



As discussed in the previous section, the two- 
dimensional Bloch band picture describing transport in 
the CPT in terms of the quasicharge leads to a 2e- 
periodic modulation of the IVC with gate voltage. This 
picture only considers Cooper pair tunneling, completely 
neglecting quasiparticle or single electron tunneling. If 
the single electron tunneling rates are large enough, they 
can disrupt this picture, causing e-periodic modulation of 
the IVC. In our experiments we find that the impedance 
of the environment can be used to suppress these quasi- 
particle tunneling rates low enough for stable 2e-periodic 
behavior to be observed. Figure shows how the pe- 
riodic modulation of the threshold voltage with gate at 
a constant current /{, = 5 pA evolves from 2e-periodic 
to e-periodic as the impedance of the environment is 
tuned. At high impedance, where the Coulomb block- 
ade is well developed, we find 2e-periodic behavior, as 
expected from the two-dimensional Bloch band picture 
(Fig. [7p). As the impedance of the environment is re- 
duced, the effect of "quasiparticle poisoning" is seen as 
a small peak in the threshold voltage at gate charges 
Qg = (2n-|-l)e (Fig.[7|3), which grows until an e-periodic 
behavior is observed (Fig. [T^V) . A similar crossover from 
2e to e-periodic behavior as the environment impedance 
is tuned has been observed by Watanabe Q and Kuo et. 
al. 113. 

To understand how quasiparticles affect the critical 
voltage of a symmetric CPT in a tunable environment it 
is instructive to determine the regions in the bias voltage 
- gate voltage plane where it is energetically favorable 
for a single quasiparticle to tunnel on or off the island. 
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Qg(e) V (^iV) 

FIG. 7: (A,B,C): Gate-induced voltage modulation of sample 
CPT 9a for increasing values of the zero bias resistance of the 
biasing SQUID arrays. The bias current is fixed at = 5 pA. 
(a,b,c): IVC of the CPT for the same values of i?o. (A, a): 
Ro « 1 Mn, (B,b): i?o « 5 Mfi, (C,c): Ro « 40 Un. 



thereby changing the island parity. As we will see, the 
shape of the corresponding stability diagrams depends on 
the impedance of the environment as well as on the ratio 
of the two characteristic energies for adding or removing 
a single quasiparticle to or from the island. These ener- 
gies are the charging energy Ec^ = e^/2Cs of the island 
and the free energy difference F — A~ kBT\ia{N^f) that 
lifts the odd parity from the even parity states 'l',^, Tg'I . 
Here, the entropy term involves iVe// = iV(0)F\/27rAT, 
the number of quasiparticle states available within tem- 
perature T just above the gap (iV(0) is the normal-state 
density of states at the Fermi level and V the volume 
of the island). At zero temperature, N^ff — and this 
energy difference is the superconducting energy gap A, 
which is paid for having one unpaired electron on the 
island. The free energy difference F vanishes at a tem- 
perature ksT* = A/ln(A^e//), generally much lower than 
the superconducting gap. 

Let us first consider the situation of a free energy differ- 
ence smaller than the island charging energy, F < Ec^ ■ 
In Fig. 8a, we plot the corresponding stability diagram 
for the case of a low impedance environment, where tun- 
neling is elastic and we can apply the global rule [23j| . 
The stability diagram contains lines which separate re- 
gions in the bias voltage - gate voltage plane with sta- 
ble parity from regions where it is energetically favorable 
for a single quasiparticle to tunnel on or off the island, 
thereby changing the parity. Within the dotted rhombi 



the transition rate Fqc from odd to even parity is zero, 
while within the solid rhombi the rate from even to odd 
parity Too vanishes. Outside the rhombi, a parity change 
corresponds to a decrease in energy and the transition 
rates Too and Too are non-zero at T = 0. We conclude 
that the odd parity states are stable near the gate charges 
Qg — [2n + l)e (Fig. 8a). At these gate charges the sys- 
tem will be "poisoned" by the excess quasiparticle which 
sits on the island in the lower energy state. When in 
the odd parity state, the effective island charge is shifted 
by e and we thus expect the gate voltage dependence of 
the critical voltage of the CPT to be shifted correspond- 
ingly (from the solid to the dashed line in Fig. 8c). In 
particular, as gate voltage is changed, the critical voltage 
should have two maxima per 2e-period (near gate charges 
Qg = (2n -|- l)e and Qg — 2ne). This is in qualitative 
agreement with the low impedance data of Fig. 7A and 
7B. 

In a high impedance environment, quasiparticle tun- 
neling is inelastic as a tunnel event is accompanied by 
ener gy exchange with the environment and the local 
rule 1231 should be applied to establish the stability di- 
agrams. As a result, the size of the rhombi doubles, 
and various stability regions start to overlap. Overlap 
means that there are bistable regions where cither parity 
is stable, depending on the history of system. However 
as can be seen in Fig. 8b, for F < Ec^, odd stabil- 
ity regions (dashed rhombi) still exist near gate charges 
Qg = (2n -I- l)e, and we expect the critical voltage to 
show two maxima per 2e period. This observation is 
in contradiction with the experimental high impedance 
data shown in Fig. 7C, where only one critical voltage 
maximum per 2e period is observed, near gate charges 
Qg = 2ne. 

We therefore consider the stability diagram for the 
case F ~ Ec^, first in a low impedance environment, 
as shown in Fig. 8d. Compared to the case F < Ec^ 
(Fig. 8a), the odd stability regions now disappear, giv- 
ing rise to regions near the gate charges Qg = {2n + l)e 
where both rates FocFoo 7^ and parity will fluctuate 
at small bias voltage. However, in a high impedance en- 
vironment for F ^ Ecs J as shown in Fig. 8e, we observe 
that the odd stability regions completely disappear, giv- 
ing way to overlapping regions of even stability only, as 
long as the applied bias voltage satisfies the inequality 
Vb < e/C — 2F/e. This is markedly different from the 
low-impedance case, where regions of fluctuating parity 
exists at any finite bias. We thus conclude that the crit- 
ical voltage is 2e-periodic as a function of gate voltage 
(Fig. 8f ) for high impedance and F ~ Ec^ , in agreement 
with the experimental findings of Fig. 7C. 

This stability analysis shows that quasiparticle poison- 
ing can occur near the gate charges Qg = (2n -I- l)e if 
conditions are not favorable. The measured, long-time 
average CPT voltage near these gate charges will depend 
on the relative weight of the odd and even states, which 
are fluctuating on a much more rapid time scale deter- 
mined by the rates FocFco. The value of these rates 
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cannot be inferred from the stability analysis alone |24j . 
For Foe ^ Feo the effective island charge will most often 
be even, hence the measured voltage of the CPT will be 
minimal near Qg = {2n + l)e (see Fig. 8f). The effects 
of poisoning will not be noticed on the long time scale 
of the experiment and one observes a stable 2e-periodic 
modulation of transport. In the opposite case Foe Foo 
near Qg = {2n + l)e, the decay of the odd state is much 
slower than that of the even state. The effective island 
charge will be predominantly odd and the measured aver- 
age voltage, instead of taking a minimum value, will reach 
a local maximum, that eventually may become equal to 
the maxima found at Qg = (2n)e (see Fig. 8c). The mea- 
surements for low-impedance presented in Fig. 7A show 
e-periodicity of the critical voltage with gate voltage, we 
conclude that the ratio Foe <C Fco. This is in agreement 
with several other experiments on CPTs with F Ecj^ in 
low impedance environments j5l l25L 12^12^ . We also note 
that the observed e-periodicity excludes that F is large 
compared Ec^ . If F > Ec^ , we would have found over- 
lapping regions of stable even parity near Qg = (2n-|- l)e, 
rather than regions of fluctuating parity and the result- 
ing dependence of critical voltage on gate voltage would 
have been 2e-periodic in the low impedance case. 

We finally comment on the effect of finite tempera- 
ture. The energy difference F depends on temperature 
T. For the islands used in this experiment, the charac- 
teristic crossover temperature T* where F ^ can be 
estimated to be approximately 250 mK. The temperature 
T appearing in the expression for F does not necessar- 
ily correspond to the cryostat temperature, and could be 
an effective temperature, describing the fluctuations in 
the environment. Previous analysis of single junctions 
IVCs biased by similar arrays 13] showed an effective 
noise temperature approximately 150 mK, which would 
result in F ^ Ec^, in agreement with our conclusions 
above. A second effect of thermal fluctuations is that 
they mix states of different parity whose energy differ- 
ence is < ksT. On the level of the stability diagrams this 
means that the lines separating stability regions of differ- 
ent parity smear into strips whose width is ~ CksT. For 
the low impedance case this has drastic effects, even at 
low bias, because it leads to broadening of the parity fluc- 
tuation regions near Qg = {2n + l)e thereby enhancing 
the e-periodic gate modulation. For high impedance the 
even parity stability region has a gap with respect to the 
bias voltage, therefore thermal fluctuations do not effect 
the 2e-periodicity as long as ksT < e^/C-2F ^ e^/2C. 
The fact that parity stability is robust against fluctua- 
tions for small bias voltage in a high impedance environ- 
ment is directly connected with the inelastic nature of 
tunneling in this limit. This inelasticity is at the origin 
of the bias voltage gap and leads to an exponential sup- 
pression of quasiparticle tunneling at low bias voltage. 
This effect is absent for low impedance as tunneling is 
elastic. 

The condition for stable 2e-periodic behavior is not 
preserved at all bias currents in the high impedance en- 
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FIG. 8: a,b,c,d: Stability diagrams for the cases F < i?Cs 
with Zcnv <C Rq and Zcnv 3> Rq (a and h), F = Ec^ with 
^cnv <S Rq and Zcnv ^ Rq (d and e). c) e-periodic gate 
voltage dependence of the critical voltage at finite bias volt- 
age (14 < e/C~2F/e) corresponding to the stability diagrams 
of a,b and d. f) 2e-periodic gate voltage dependence of the 
critical voltage at finite bias voltage (14 < e/C — 2F/e) cor- 
responding to the case e). 



vironment. Figure shows the gate voltage modulation 
for high impedance {Rq « 20 Mf]) of the CPT voltage 
for three different values of the bias current. ^From low 
current levels up to the back-bending region of the IVC, 
the voltage modulation is stable 2e-periodic (Fig. 13^). As 
the current is increased toward the Zener current, which 
for this sample is measured at Iz ~ 75 pA, the rate of 
single electron transitions increases, as explained in the 
previous section. Zener tunneling causes an increase in 
the rate Foo near Qg = (2n-|-l)e, and the effect of poison- 
ing can be observed (Fig. 03). At even higher currents, 
the quasiparticle tunneling rate is so large that only e- 
periodic modulation can be seen (Fig. |5p). 



V. CONCLUSIONS 

The effect of the electrodynamic environment on single 
charge tunneling in the Cooper Pair Transistor (CPT) 
has been studied experimentally by biasing and mea- 
suring the CPT through one-dimensional SQUID arrays. 
The impedance of these arrays can be tuned with an ex- 
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arrays is increased, the CPT develops a Coulomb block- 
ade of Cooper pair tunneling with a well defined back- 
bending current- voltage characteristic (IVC), indicative 
of the coherent transfer of single Cooper pairs. The gen- 
eral shape of the measured IVC is that predicted by a 
theory based on quasicharge dynamics. We find that the 
observed threshold voltage depends on the ratio between 
the Josephson energy and the charging energy in a way 
that is in qualitative agreement with the Bloch band the- 
ory for the CPT. As the impedance of the environment is 
increased, we observe a transition from e periodic to 2e- 
periodic modulation of the IVC, demonstrating that the 
high impedance environment restores the parity effect 
in the CPT. Enhanced parity effect, or reduced quasi- 
particle poisoning, is explained by the inelastic nature 
of single-charge tunneling in a high impedance environ- 
ment, where a charging energy gap makes odd parity 
states unstable near gate charges Qg = (2n + l)e. 
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